Stellar Populations - Planning for the Next Decade 
Proceedings IAU Symposium No. xxx, 2009 
G. Bruzual, & S. Chariot, eds. 



© 2009 International Astronomical Union 
DOI: 00.0000/X000000000000000X 



Nebular abundances in galaxies: 
Beware of biases 

Grazyna Stasinska 1 

1 LUTH, Observatoire de Paris, CNRS, Universite Paris Diderot; Place Jules Janssen 92190 

Meudon, France 
email: grazyna. stasinska@obspm.fr 

' Abstract. The derivation of nebular abundances in galaxies using strong line methods is simple 
and quick. Various indices have been designed and calibrated for this purpose, and they are 
widely used. However, abundances derived with such methods may be significantly biased, if 
the objects under study have different structural properties (hardness of the ionizing radiation 
field, morphology of the nebulae) than those used to calibrate the methods. Special caution is 
required when comparing the metallicities of different samples, like, for example, blue compact 
galaxies and other emission line dwarf galaxies, or samples at different redshifts. 
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U: 1. Introduction 

. Why talk about nebular abundances in a symposium devoted to stellar populations in 
galaxies? The ultimate goal of stellar populations studies is to understand the evolution 
Q . of galaxies. In principle, the determination of the metallicities of the stellar populations 
gives an evolutionary view of the chemical enrichment of galaxies. Such studies are just 
beginning (see e.g. Vale Asari et al. 2009 and contribution to this symposium). For the 
i moment, however, the overwhelming majority of studies on the chemical composition of 

galaxies rely on global stellar metallicity indices (for early- type galaxies), and on nebular 
abundances (for star-forming galaxies) derived from emission-line studies. Note that the 
abundances derived by one or the other method do not have the same meaning: Stellar 
indices provide an integrated (and weighted) metallicity of the different generations of 
stars while the chemical composition of nebulae is the result of mixing with the interstellar 
gas of the metals ejected by all the present and past stellar generations and is the same 
. as that of the presently forming stars. 
t-H ' In this communication, we would like to draw attention to widely overlooked biases in 
ON . the determination of nebular abundances. We will restrict to the so-called "strong line 
methods" for abundance determination (see e.g. Stasinska 2004 for a general introduction 
. on nebular abundance determinations). After presenting a number of abundance calibra- 
tors, We will show how their use may, in certain cases, lead to erroneous statements. 
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2. A summary of strong line methods and their recent calibrations 

Strong line methods are statistical in nature, and take advantage of the fact that, in 
giant HII regions, the metallicity appears to be linked to the mean effective temperature 
of the stellar radiation field and to the ionization parameter. They are presently being 
applied to integrated spectra of galaxies (e.g. Kewley & Ellison 2008), which poses the 
yet unsolved problem of how the galactic abundance gradients affect the results. Here, we 
will ignore this problem, and show that abundances obtained by strong line methods can 
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be biased even in the case of chemically homogeneous metal-poor galaxies and, therefore, 
lead to incorrect interpretations. 

Strong line methods, being statistical, have to be calibrated. This can be done either 
using a grid of photoionization models or a sample of objects for which the oxygen 
abundances have been determined independently by direct methods. 

Here is a list of abundance calibrators which are routinely used for nebular abundance 
studies. 

• a 23 , i-e. ([On] A3727 + [Om] A5007)/H/?, introduced by Pagel et al. (1979) 

• b 3 N 2 , i.e. [Om] A5007/[Nn] A6584, proposed by Alloin ct al. (1979) 

• c S 23 , i.e. ([Sn] A6716,A6731 + [Sm] A9069)/H/3, proposed by Vilchcz & Esteban 
(1996) 

• d N 2 , i.e. [Nil] A6584/Ha, first used by Storchi Bergmann et al. (1994) 

• e S3O3, i.e. [S in] A9069/[Om] A5007, proposed by Stasiriska (2006) 

• f Ar 3 3 , i.e. [Arm] A7135/[Om] A5007, proposed by Stasiriska (2006) 

• g Ne 3 2 , i.e. [Nem] A3869/[Om] A5007, proposed by Nagao et al. (2006) 

• h O3O2, i.e. [Om] A5007/[On] A3727, considered by Nagao et al. (2006) 

As emphasized by Stasiriska (2008), not all of those indices have a direct link with 
the metallicity (or rather with the oxygen abundance, which is taken as a proxy for 
metallicity in nebular studies). It is the empirical relation between metallicity and struc- 
tural properties of HII regions that allows them to be used as metallicity indicators. For 
example, the hardness of the ionizing radiation in giant HII regions is correlated with 
metallicity, due to opacity effects in stellar interiors atmospheres that soften the ionizing 
radiation field (Mc Gaugh 1991). The ionization parameter is correlated with metallicity 
due to stellar wind pressure that shapes the morphology of the emitting gas (Dopita et 
al 2006). Finally, N/O is related to O/H due to secondary production of nitrogen. 

Those indices (and others) have been calibrated many times. Here, we will use the 
calibration by Nagao et al. (2006) for O23, O3N2, N2, Ne302 and O3O2, Pcrcz-Montero 
& Diaz (2005) for S 23 and Stasihska (2006) for S3O3 and Ar 3 3 . 



3. Evidence for biases with strong line methods 

The following is motivated by the strange result of Sanchez Almeida et al. (2008) that 
BCDs (blue compact dwarf galaxies) in the Sloan Digital Sky Survey (SDSS) have, on 
average, lower oxygen abundances than QBCDs (quiet blue compact dwarf galaxies). 
This claim is based on abundances derived using the N2 index. 

That oxygen abundances derived from methods involving nitrogen lines may be biased 
has been noted many times. We will show below that the problem is more serious than 
that and affects any statistical method. We consider all the emission line galaxies from the 
Sloan Digital Sky Survey (data release 6, Adelman-McCarthy et al. 2008) for which the 
weak [O in] A4363 line has been measured with better than 30% accuracy, and divide this 
sample in two groups. One which has H/3 equivalent widths, Wup, larger than 100A, the 
other with smaller values of Wu/3- The first group bears similarities with the BCD class 
of Sanchez Almeida ct al. (2008) (newborn galaxies), while the second group contains 
objects more like their QBCDs (aging galaxies), since low values of Wa/3 tend to indicate 
more evolved populations of the ionizing stars. For all those objects, we can obtain the 
oxygen abundance with the direct, temperature-based method, following the recipes given 
in Izotov et al. (2006). We can also compute the oxygen abundances using strong line 
indices listed above. In the various panels of Fig. 1, we plot the values of those indices 
as a function of O/H as derived from the direct method. BCDs are represented by black 
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Figure 1. The values of some strong-line ratios used as abundance indicators, as a function of 
the oxygen abundance. The curves represent various calibrations (as indicated in the text). The 
dots represent galaxies from the SDSS for which O/H was derived with a temperature-based 
method. Black dots represent galaxies with Wh/3 > 100A. Cyan dots represent galaxies with 
Wh/3 < 100A. Clearly, strong line methods tend to underestimate the values of O/H for black 
points and overestimate them for cyan points. 



points, QBCDs by cyan (light grey) points. The calibrating functions presented in the 
references at the end of Sect. 2 are shown as continuous lines. 

The result is eloquent. QBCDs are vertically displaced from BCDs in all the diagrams. 
There is no indication of a difference in O /H values (derived with direct methods) between 
the two groups. On the other hand, if one were to use the strong line indices, Fig. 1 clearly 
shows that the O /H values would tend to be underestimated for BCDs and overestimated 
for QBCDs. 

Thus, the application of inadequate abundance calibrators to two different samples of 
objects leads to the erroneous conclusion the the two samples have different mctallicitics. 
We believe that it is exactly what is happening in the case of the analysis by Sanchez 
Almeida et al. (2008, 2009). 

Why is it so? QBCDs, having older ionizing stellar populations than BCDs, are ionized 
by a softer radiation field. In addition, their HII regions are presumably less compact 
and more filamentary because they have suffered the mechanical action of stellar winds 
for a longer time than BCDs. As a consequence, their ionization parameter is smaller. 
Both reasons contribute to lowering the mean excitation of QBCD HII regions with 
respect to BCD ones. All the calibrators used here (except O23, in panel a) are strongly 
dependent on the ionization state of the gas. Therefore, they cannot serve to compare 
the metallicities of two samples which systematically differ in their ionization state! We 
note that, unfortunately, all the calibrators are biased in the same direction, therefore 
using several calibrators will not solve the problem. 

The O23 calibrator does not depend strongly on the ionization state, but it does depend 
on the hardness of the heating radiation. Therefore, if is not fully secure either (and, as 
is known, it is double valued which makes its use complicated). 
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4. Conclusion 

Strong line methods to derive nebular abundances are very easy to apply but they 
are prone to systematic errors. In principle, they should be used only for objects whose 
HII regions have the same structural proprties as those of the calibrating samples. This 
recommendation is not easy to follow, but at least one should be aware that using the 
same calibration for different samples may produces important biases. 

In particular, claims on differences in oxygen abundance 

• between samples of galaxies with different chemical evolution histories 

• between samples of galaxies with different star formation histories 

• between samples of galaxies at different redshifts (observational selection effects may 
play a role) 

should be taken with, at a minimum, a grain of salt. 
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